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Geminga is the second brightest persistent source in the GeV 7-ray sky. Discovered in 1975 
by SAS-2 mission, it was identified as a pulsar only in the 90s, when ROSAT detected the 237 
ms X-ray periodicity, that was later also found by EGRET in 7 rays. Even though Geminga has 
been one of the most intensively studied isolated neutron star during the last 30 years, its interest 
remains intact especially at 7-ray energies, where instruments like the Large Area Telescope (LAT) 
aboard the Fermi mission will provide an unprecedented view of this pulsars. We will report on the 
preliminary results obtained on the analysis of the first year of observations. We have been able 
to do precise timing of Geminga using solely 7 rays, producing a timing solution and allowing a 
deep study of the evolution of the light curve with energy. We have also measured and studied the 
high-energy cutoff in the phase-averaged spectrum and produced a detailed study of the spectral 
evolution with phase. 



I. INTRODUCTION 

In the last 30 years Geminga was for a long time the 
first and only radio-quiet 7-ray pulsar detected in the 
sky. After Fermi launch it is no more alone, for other 
7-ray pulsars have been discovered without a radio or 
X-ray counterpart so far. Since its discovery Geminga 
has been widely and deeply studied, but nowdays its 
importance is increased since it is the harbinger of a 
population of radio-quiet neutron stars. Geminga was 
discovered in 1975 as a 7-ray source by SAS-2 [lj] 
but only in 1992 it was identified as a pulsar object 
thanks to the ROSAT discovery of a periodic X-rays 
emission from this source [16| . This has allowed a 
successful search for periodicity in high-energy 7 rays 
with EGRET Q . Recent parallax and proper motion 
measurements confirm the earlier results, yielding a 
distance of 250_^° pc and a proper motion of 178.2 
± 0.4 mas/yr (13j ]. Geminga has a period of 237 ms 



and a period derivative of 1.1 
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these 



parameters imply that it is an old 7-ray with a char- 
acteristic age of 3 x 10 5 yr and spin-down luminosity 
E = 3.26 x 10 34 erg s _1 . The light curve detected by 
EGRET showed two peaks separated by almost half 
rotation [llj . The detected spectrum was compat- 
ible with a simple power law with a photon index of 
1.42 even if there was evidence of a cutoff at ~ 2 GeV, 
but the limited EGRET statistics did not allow a mea- 
surement. Deep X-ray observations allowed XMM- 
Newton and Chandra to map the neutron star surface 
as it rotates, bringing into view different regions con- 
tributing different spectral components [1, [HI EH as 



well as an arcmin-scale bow-shock feature trailing the 
pulsar's motion @, EH- Such a non-thermal diffuse 
X-ray emission surrounding the pulsar is produced by 
high-energy (~ 10 14 eV) electrons escaping from the 
pulsar magnetosphere and interacting with the inter- 
stellar magnetic field. 

Fermi is an excellent instrument for high energy 7-ray 
pulsar studies, thanks to its outstanding timing capa- 
bility, really good pointing and energy resolution. In 
this proceeding we present the analysis of the Geminga 
pulsar based on the statistics collected during the first 
year of operations of the Fermi mission. More results 
on the Geminga observations by Fermi can be found 
ind. 



II. GAMMA-RAY OBSERVATIONS 

Fermi-LAT is a pair conversion telescope sensitive 
to 7 rays from 20MeV to >300GeV. It is composed 
by a silicon tracker, a Csl hodoscopic calorimeter and 
an Anticoincidence detector to reject the charged par- 
ticles background The LAT has a large effective 
area (~8000 cm 2 on-axis) and thanks to its large field 
of view (~ 2.4 sr), it covers the entire sky every 2 
orbits (~ 3 h). These characteristic lead to an in- 
crement of a factor ~ 30 in sensitivity with respect 
to its predecessor EGRET. The data used are of par- 
ticular importance since cover the full first year of 
the Fermi mission. The collected data range from 
Launch and Early Operations (L&EO) (covering ~ 
two months after 2008 June 25, when the LAT was op- 
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erated in pointing and scanning mode for tuning pur- 
poses) to the first year of nominal operations (namely 
2009 June 15). Only for the spectral analysis we se- 
lected a subsample of data, covering observations in 
scanning mode under nominal configuration from 2008 
August 4 to 2009 June 15. We excluded observations 
when Geminga was viewed at zenith angles > 105° to 
the instrument axis where Earth's albedo 7 rays give 
excessive background contamination. We also have 
taken care of excluding time intervals when the Re- 
gion Of Interest (ROI) intersects the Earth's albedo 
region. 

III. RESULTS 
A. Pulsed profile 

The timing solution used to build the light curve 
has been obtained using only the LAT data, thanks 
to its excellent timing accuracy (better than 1/xs). We 
have used TEMP02 [TtJ to generate a pulse profile 
and then we have fitted the results with the timing 
solution model: the residuals have an RMS of 251/is, 
as shown in Figure [TJ Like all the pair conversion 
telescopes also the Fermi-LAT has a Point Spread 
Function (PSF) that strongly depend with energy. 
We have chosen to select photons coming from 
Geminga with an angle $ <max[1.6°-3° logio(EGev), 
1.3°], that maximize the pulsed signal over the 
background. The light curve is shown in Figure O 
It has been plotted using a variable-width phase 
bins, each one containing 400 events, to enhance 
the fine structure. The photon flux in each phase 
interval, i.e. the counts divided by phase width, 
have thus a la Poisson statistical error of 5%. The 
dashed line represents the contribution of the diffuse 
background, as estimated by selecting photons in the 
phase interval <j) = 0.9-1.0 from an annular region 
within 2° and 3° from the source. The total number 
of pulsed photons is 61219 ± 247 and the backgorund 
photons are 9821 ± 99. Both the peaks have been 
fitted with an asymmetric Lorentzian function and 
we obtain the maximum for peak 1 at <p = 0.141 
± 0.002 and for peak 2 at <p = 0.638 ± 0.003: the 
separation is 0.497±0.004. We have also built light 
curves in different energy bands (Figure [3]) and fitted 
the peaks in each of them. The first trend that can be 
seen is the first peak is decreasing at higher energies, 
confirming a decreasing trend in the ratio P1/P2 as 
observed in Crab, Vela and PSR B1951+32 7-ray 
pulsars by EGRET [23[ and now confirmed for many 
pulsars, including the Vela [l| and the Crab pulsars 
[3| by Fermi LAT. We observed pulsed photons up to 
- 18 GeV. 

Moreover, both peaks narrow with increasing energy. 
The decreasing trend in pulse width of PI and P2 is 
nearly identical. The emission of the "first interpeak" 



between PI and P2 is significantly detected up to 
10 GeV, emission in the "second interpeak" region 
(between 0.9 and 1.0), not detected before, is clearly 
present at low energies but vanishes above ~ 2 GeV. 



B. Energy Spectrum and phase-resolved 
spectroscopy 

The spectral analysis is one of the key and most 
powerful analysis for pulsars. To perform the spec- 
tral analysis we have used the maximum-likelihood 
estimator gtlike, that is part of the standard Fermi 
Science Tools provided by the FSSC[25]. The fit was 
performed using a region of the sky 15° around the 
pulsar position between 0.1 and 100 GeV. To prop- 
erly fit the source we have taken into account all the 
nearby source in a 20° around Geminga, the galactic 
diffuse emission based on GALPROP model [22J and 
the isotropic residuals emission based on a model de- 
rived from a fit of LAT data at high latitude. 
The function that we have used to fit the spectral 
emission is a power law with exponential cutoff in the 
form: 

^ = N E~ r exp (--|) cm-^GeV- 1 (1) 

The fit gives the following spectral parameters No — 
(1.19 ± 0.01 ± 0.08) x 10~ 6 cm" 2 s" 1 GeV" 1 , T = 
(1.30 ± 0.01 ± 0.05) and E = (2.47 ± 0.04 ± 0.19) 
GeV. First error is statistical and the second is sys- 
tematic. The integral flux above 0.1 GeV is (4.15 ± 
0.02 ± 0.39) x 10~ 6 cm" 2 s" 1 and the correspond- 
ing energy flux is (2.57 ± 0.01 ± 0.32) x 10~ 3 MeV 
cm~ 2 s _1 . We have also fitted the data with a simple 
power law and with a power law with a super expo- 
nential cutoff with the super exponent fixed at 2, and 
we have rejected both at more then 5a level. 
Thanks to the huge number of photons collected we 
were able also to perform phase-resolved spectroscopy. 
We divided the period in 36 phase bins, each contain- 
ing 2000 counts, according to the energy-dependent 
cut presented in Sec. IIII Al We have found that the 
emission is always well described by a power law with 
exponential cutoff, even if with different spectral pa- 
rameters. This means that the emission in the two 
peaks is different, as expected from the study of the 
light curve with energies, and that the pulsar emis- 
sion is clearly visible in both the interpulse and in the 
off-pulse region. The phase evolution of the photon 
index and the cutoff energy (Figure 0|) show a clear 
variation across the whole rotation. One of the main 
features is that the "second interpeak" in the range <j> 
~ 0.9-1.0 shows an exponential cutoff, indicating that 
the pulsar emission extends over the whole rotation. 
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FIG. 1: Timing solution of Geminga obtained solely from 7 rays. (The Y-scale is chosen to highlight the timing noise 
residuals, that corresponds to only 0.001 in phase) 



IV. DISCUSSION 

One of the key feature of Geminga is the absence 
of radio emission, that clearly favor the outer mag- 
netospheric emission models. The current evidence 
against low-altitude emission in 7-ray pulsars [lj can 
also be supplemented by constraints on a separate 
physical origin. The maximum energy of the pulsed 
photons observed must lie below the 7-B pair produc- 
tion mechanism threshold, providing a lower bound 
to the altitude of the 7-ray emission r m i n . According 
to [||, the highest energy of pulsed photons observed 
from Geminga (E~18 GeV), points toward a limit of 
"Train > 2.7 stellar radii, a value clearly precluding 
emission very near the stellar surface, adding to the 
advocacy for a slot gap or outer gap acceleration lo- 
cale for the emission in this pulsar. 
Following the atlas of 7-ray light curves compiled by 
we can use Gemingas light curve to estimate, for 
each model, the stars emission parameters, namely the 
Earth viewing angle Qe with respect to the neutron 
star spin axis, and the inclination angle a between the 
stars magnetic and rotation axes. For Slot Gap/Two 
Pole Caustic model [12l . 1201 ] we obtain a in the range 
of 30 - 80, 90 degrees and C in the range 90, 55 - 80 
degrees, while for Outer Gap [2l| models we obtain a 
in the range of 10 — 25 degrees and ( ~ 85 degrees. 
The total luminosity radiated by the pulsar is then 
given by L 1 = 4ir fnF b s D 2 where F f, s is the ob- 



served phase-averaged energy flux over 0.1 GeV and 
D = 2507^120 P c * s tne P u l sar distance [l3[. The beam- 
ing factor /n can be derived from [24| and points to a 
value of 0.7 - 0.9, 0.6 - 0.8 for Slot Gap and 0.1 - 0.15 
for Outer Gap. The estimated averaged luminosity 
is then L 7 = 3.1x10 34 /q erg s _1 , yielding a 7-ray ef- 
ficiency 77 7 = ^ = 0.15/n (d/100pc) 2 . Altough the 
uncertainty in the distance limit the conclusions, the 
phase-resolved spectral variation observed can be pre- 
dicted by caustics models. 



V. SUMMARY 

We presented the analysis of Geminga based on data 
collected during the first year of operations of Fermi. 
We were able to obtain a timing solution based only 
on 7 rays and to produce a high-detailed light curve, 
that shows clear energy-dependent variation of the 
profile. Although the phase-averaged spectrum is con- 
sistent with a power law with exponential cut-off, the 
phase-resolved analysis showed a much richer picture 
of different spectral components intervening at differ- 
ent rotational phases. The phase-resolved analysis has 
also allowed the detection of the "second interpeak" 
emission indicating a pulsar emission extending over 
all phases. Our results favor the outer magnetospheric 
origin for the 7-ray emission and future improvements 
in estimating the distance of Geminga will help to bet- 
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FIG. 2: Light curve of Geminga above 0.1 GeV. The profile has been plotted in fixed- width phase bins of 400 counts 
each to highlight fine details. The insets show details of the two peaks and the interpulse 



ter strenghten the conclusions and constraining outer 
magnetospheric models. 
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FIG. 4: Phase evolution of the photon index (top) and energy cut-off (bottom) above 0.1 GeV as the function of the 
pulse phase, divided in phase bins each containing 2000 photons. The dashed histogram represents the Fermi-LAT light 
curve above 0.1 GeV in variable- width phase bins of 2000 photons/bin. 
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